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Abstract 

We review our recent results on short time approxima- 
tions, with emphasis on applications for which the system- 
environment interactions involve a general non-Hermitian 
system operator A, and its conjugate, . We evaluate the 
onset of decoherence at low temperatures in open quan- 
tum systems. The developed approach is complementary 
to Markovian approximations and appropriate for evalua- 
tion of quantum computing schemes. Example of a spin 
system coupled to a bosonic heat bath via A oc cj_ is dis- 
cussed. 



1 Introduction 

The coupling of a quantum system to environmental de- 
grees of freedom induces decoherence, destroying quan- 
tum superposition and reducing pure states to mixed 
states. Understanding decoherence is important for quan- 
tum control and computing and, generally, for obtaining a 
description to the evolution of the system's reduced den- 
sity matrix with the environmental modes traced over. 
Since this can not be done exactly in most cases, differ- 
ent approximation techniques which are valid for differ- 
ent time scales were developed. For short times, appropri- 
ate for quantum computing gate functions and, generally, 
for controlled quantum dynamics, approximation schemes 
for the density matrix have been suggested recently [1^-]. 
The present survey introduces a new 1 1 1 rather general 
short-time approximation which applies for models with 
system-bath interactions involving a general system oper- 
ator. It thus extends the previously known approach [3,4] 
developed for couplings involving a single Hermitian sys- 
tem operator. 

We assume that the Hamiltonian of the open quantum 
system is 

H = H S + H B + H I , (1) 

where Hs describes the system which is coupled to a fluc- 
tuating dynamical reservoir (the bath). Typically, the bath 



is modeled by the harmonic modes, as reviewed in 1 5 1, 

H B =J2"kblh- (2) 

k 

Here are the annihilation operators of the bath modes, 
and we use the convention % = 1. We assume that the in- 
teraction with the bath involves the system operator A that 
couples linearly to the bath modes, as reviewed in |6|, 

k k 

with the interaction constants g k . 

Let R(t) denote the overall density matrix. It is as- 
sumed [3,7] that at time t — the system and bath 
are not entangled, and the bath modes are thermalized 

(J3 = 1/kT): 

R(0) = p(0)l[9 k , (4) 

k 

Q k = {l-e-^ k )e-^ h l^. (5) 

We point out that while the quantum system S, de- 
scribed by the reduced density matrix p(t) = Tr^i? (t), 
is small, typically two-state (qubit) or several-qubit, the 
bath has many degrees of freedom. The combined effect 
of the bath modes on the system can be large even if each 
of them is influenced little by the system. This has been 
the basis for the arguments for the harmonic approxima- 
tion for the bath modes and the linearity of the interac- 
tion, as well as for the Markovian approximations [6,7] 
that assume that the bath modes are "reset" to the thermal 
state by the "rest of the universe" on time scales shorter 
than any dynamical time of the system interacting with 
the bath. 

The frequencies of the oscillators of the bath are usually 
assumed to be distributed from zero to some cutoff value 
cu c . The bath modes with the frequencies close to the en- 
ergy gaps of the system, AEij — Ei — Ej, contribute to 
the "resonant" thermalization and decoherence processes. 
Within the Markovian schemes, the diagonal elements of 
the reduced density matrix of the system, approach the 
thermal values oc eT Ei l kT for large times exponentially, 
on time scale T%, The off-diagonal elements vanish, which 



represents decoherence, on time scale Ta, which, for res- 
onant processes, is given by T2 ~ 27i. However, gen- 
erally decoherence is expected to be faster than thermal- 
ization because, in addition to resonant processes, it can 
involve virtual processes that do not conserve energy. It 
has been argued that this additional "pure" decoherence is 
dominated by the bath modes with near-zero frequencies 
[3,6,7]. At low temperatures, this "pure decoherence" is 
expected QD to make T 2 -C T\. 

Since the resetting of these low-frequency modes to the 
thermal state occurs on time scales fi/kT — (3, the Marko- 
vian approach cannot be used at low temperatures [3,7]. 
For quantum computing in semiconductor-heterostructure 
architectures [8-15], temperatures as low as few 10 mK 
are needed. This brings the thermal time scale to (3 ~ 
10~ 9 sec, which is close to the single-qubit control times 
10~ 11 -10 -7 sec. Alternatives to the Markovian approxi- 
mation have been suggested [1^4,16-21]. 

2 Short-Time Approximation 

In applications in quantum computing, calculations with 
only a single qubit or few qubits are necessary for evalu- 
ation of the local "noise," to use the criteria for quantum 
error correction [22-27]. For example, the system Hamil- 
tonian is frequently taken proportional to the Pauli matrix 
a z . The interaction operator A can be proportional to a x , 
which is Hermitian. Such cases are covered by the short- 
time approximation developed earlier [3,4]. However, one 
can also consider models with A oc ct_ . Similarly, models 
with non-Hermitian A are encountered in Quantum Optics 
1281 . In this section, we develop our short time approxi- 
mation scheme. Results for a spin-boson type model are 
given in the next section. 

We derive a general expression for the time evolution 
operator of the system within the short time ap- 

proximation. The overall density matrix, assuming time- 
independent Hamiltonian over the quantum-computation 
gate function time intervals [8-15], evolves according to 

R(t) = U(t)R(0)[U(t)?, (6) 

U(t) = e -KHs+H B+Hl )t_ (?) 

The general idea of our approach is to break the exponen- 
tial operator in Q into products of simpler exponentials. 
This involves an approximation, but allows us to replace 
system operators by their eigenvalues, when spectral rep- 
resentations are used, and then calculate the trace of R(t) 
over the bath modes, obtaining explicit expressions for the 
elements of the reduced density matrix of the system. For 
Hermitian coupling operators, A T = A, our approach re- 
duces to known results [3,4]. 



In order to define "short time," we consider dimension- 
less combinations involving the time variable t. There 
are several time scales in the problem. These include 
the inverse of the cutoff frequency of the bath modes, 
1 /oj c , the thermal time (3=1/ kT, and the internal char- 
acteristic times of the system 1/A.Ey. Also, there are 
time scales associated with the system-bath interaction- 
generated thermalization and decoherence, T\ 3. The 
shortest time scale at low temperatures (when (3 is large) 
is typically l/u/ c . The most straightforward expansion in 
t yields a series in powers of uj c t. The aim of developing 
more sophisticated short-time approximations [1,3,4] has 
been to preserve unitarity and obtain expressions approxi- 
mately valid up to intermediate times, of order of the sys- 
tem and interaction-generated time scales. The applicabil- 
ity for intermediate times can only be argued for heuristi- 
cally in most cases, and checked by model calculations. 

We split the exponential evolution operator into terms 
that do not have any noncommuting system operators in 
them. This requires an approximation. For short times, 
we start by using the factorization 

e -i(H s +H B +H,)t+0(t 3 ) 
= e -$Hst e -i(Hi+H B )t e -§H 8 t (g) 

where we have neglected terms of the third and higher or- 
ders in t, in the exponent. The middle exponential in {8J, 

g = e -i(H B + H,)t _ e - l (H B +AG t +A t G)t ^ 

where 

G = J29pk, (10) 

k 

still involves noncommuting terms as long as A is non- 
Hermitian. In terms of the Hermitian operators 

L = i(A + At), (11) 

M = l - (A - At) , (12) 

we have 

AG t + A t G = L(G + G t )+iM(G-G t ). (13) 

We then carry out two additional short-time factorizations 
within the same quadratic-in-i (in the exponent) order of 
approximation, 

3 = e h[M(G-G^)-iH B ]t e iH B t (U) 
x e -i[H B +L(G+G'<)}t e i [ H B t e ±[M(G-G' t )-iH B }t 

This factorization is chosen in such a way that 5 remains 
unitary, and for M = or L = the expression is identi- 
cal to that used for the Hermitian case [3,4]. The evolution 
operator then takes the form 



with S from dl4t . which is an approximation in terms of a 
product of several unitary operators. 

It has been recognized [1,3,4] that approximations of 
this sort, which are not perturbative in powers of Hi, are 
superior to the straightforward expansion in powers of 
uj c t. Specifically, in (|8jl, we notice that H$ is factored out 
in such a way that Hb, which commutes with Hs, drops 
out of many commutators that enter the higher-order cor- 
rection terms. This suggests that a redefinition of the ener- 
gies of the modes of Hb should have only a limited effect 
on the corrections and serves as a heuristic argument for 
the approximation being valid beyond the shortest time 
scale 1/cj c , up to intermediate time scales. 

Our goal is to approximate the reduced density matrix 
of the system. We consider its energy-basis matrix ele- 
ments, 

Pmn (i) = Tr B (m| UR (0) C/ t \n) , (16) 

where 

H s \n) = E n \n) . (17) 

We next use the factorization (|8j, dl4t to systematically 
replace system operators by c-numbers, by inserting de- 
compositions of the unit operator in the bases defined by 
Hs, L, and M. First, we collect the expressions ©, dl4t . 
dl5> . dl7> . and use two energy-basis decompositions of 
unity to get 

Pmn (t) = ]Tei^+^«*)Wo) 

x Tr s [HH|p) Jj0 fc <g|S+|n>].(18) 

k 

We define the eigenstates of L and M, 

L\\) = \ | A) , M \n) = jj, \fi) . (19) 

The operators S and St introduce exponentials in (1181 that 
contain L or M in the power. By appropriately insert- 
ing 2~2 l A ) (M or 2~2 Im) (tA between these exponentials, 

we can convert all the remaining system operators to c- 
numbers. 

Now the trace in Jl 8i can be evaluated, by using opera- 
tor identities for bosonic operators 1 28 1 and the coherent- 
states technique. We obtain our final result for the density 
matrix evolution 1 1 1, 

p,q fijXj 

x (m\ Mi) (Mil h) (Ai |Ma) (M2I p) 

x (q\ fx 3 ) (/ia| A 2 ) (A 2 \m) <A* 4 ] n) , (20) 

where the first sum over p and q is over the energy eigen- 
states of the system; the second sum is over Ai,A 2 and 



Ml, ••• , M4> which label the eigenstates of the operators L 
and M, respectively. The power in the exponential is 

V = B 2 (t)(\ 2 _+ f /_^)+B 2 (t/2)(^-p'_) 2 

- F (t) (m" - M-) A - - iC (*) A - A + - iC (t/2) 
x (mV+ + M-M+) + i5 (t) (A_m+ - A+M_) 

- id(»MlM+- (21) 
Here we introduced the variables 

M± = Mi±M4, (22) 
M± = M2±M3, (23) 

and 

A ± =A!±A 2 . (24) 
and the spectral sums over the bath modes, 

5^ (t) = 2 > — — sm ~ o - coth ~W > ( 25 ) 

C (t) = V M- (wjtt - sin w fc t) ; (26) 
k 

these functions are well known [29,30]. The result also 
involves the new spectral functions 

S(t) = -2£^sin 2 ^, (27) 

FW = 4 y^sin 2 ^sin^coth^. (28) 
k Wfe 4 2 2 

Furthermore, we defined 

d (t) = 2C (t/2) -C(t). (29) 

3 Discussion and Application 

In most applications evaluation of decoherence will re- 
quire short-time expressions for the reduced density ma- 
trix of a single qubit. Few- and multi-qubit systems will 
have to be treated by utilizing additive quantities [31-33], 
accounting for quantum error correction (requiring mea- 
surement), etc. For a two-state system — a qubit — the 
summation in i20\ involves 2 8 = 64 terms, each a product 
of several factors calculation of which is straightforward. 
We carry out the calculation for an illustrative example. 
We consider the model 1 34 1 defined by 

H = A o z + "Ah + E {^A + 9* k <r+h) , 
k k 

(30) 



where A > is a constant, <r± = \{a x ± ia y ) and 
er 2 are the Pauli matrices, bj, and 6 fe are the bosonic cre- 
ation and annihilation operators, and g). are the coupling 
constants. Physically this model may describe, for exam- 
ple, a qubit interacting with a bath of phonons, or a two- 
level molecule in an electromagnetic field. In the latter 
case, this is a variant of the multi-mode Jaynes-Cummings 
model [28,35]. Certain spectral properties of this model, 
the field-theoretic counterpart of which is known as the 
Lee field theory, are known analytically, e.g., |36|. How- 
ever, the trace over the bosonic modes, to obtain the re- 
duced density matrix for the spin, has not been calculated 
exactly. 

For the model POt we have A = cr_ and At = <j + , 
so that L = a x /2 and M = <J y /2. We have | Ai 5 s) = 
(|t) ± II)) / V2, with the eigenvalues Ai, 2 = ±1/2, and 
1/^1,2) = (IT) ± i ID) /\/2, with the eigenvalues ^ ii2 = 
±1/2. For the initial state, let us assume that the spin at 
t = is in the excited state |D (f |, so that the initial den- 
sity matrix has the form 



P(0) 



1 




(31) 



Our calculation yielded the following results for the den- 
sity matrix elements: pi2(t) = and 

4pu(i) = 2+ e- 2B2 ( t )+e- 4B2 (*)cosh(2F) 

+ 2e~ 2B2 ^ sinh {B{) cos (S) + 2e~ s2( ^) 

x cos (d) sin (S^-^-^) [e* Cl 

x sinh {-iS + F)+ e- lCl sinh {-iS - F)] 

where C\ was defined in \29l and 

B Y (t) = 2B 2 {t/2) - B 2 (t) . (32) 

Where not explicitly shown, the argument of all the spec- 
tral functions entering (I32> is t. 

In order to obtain irreversible behavior and evaluate 
a measure of decoherence, we consider the continuum 
limit of infinitely many bath modes. We introduce the 
density of the bosonic bath states T> (uS), incorporating 
a large-frequency cutoff uj c , and replace the summations 
in d25i — (I28t by integrations over li [5,29,37,38]. For in- 
stance, d25i takes the form, 



uj 2 2 2 



B 2 (t) = J ' <L> 



We will use the standard Ohmic-dissipation expres- 
sion, with an exponential cutoff, for our illustrative calcu- 
lation, 



where fl is a constant. 

Our results for the density matrix elements depend on 
the dimensionless variable oj c t, as well as on the dimen- 
sionless parameters £1 and u> c f3 (= Tiuo c /kT, where we 
remind the reader that %, set to 1, must be restored in 
the final results). Interestingly, the results do not de- 
pend explicitly on the energy gap parameter A, see OOt . 
This illustrates the point that short-time approximations 
do not capture the "resonant" relaxation processes, but 
rather only account for "virtual" relaxation/decoherence 
processes dominated by the low-frequency bath modes. 
However, the short-time approximations of the type con- 
sidered here are meaningful only for systems with well- 
defined separation of the resonant vs. virtual decoherence 
processes, i.e., for h/A 1/oj c . For such systems, 
h/A= l/A defines one of the "intermediate" time scales 
beyond which the approximation cannot be trusted. 

As an example, we calculated a measure of deviation of 
a qubit from a pure state in terms of the "linear entropy" 
[31,33,39], 



,s(<) = l-Tr[p 2 (<)] 



(35) 



(34) 



Figure 1 schematically illustrates the behavior of s(i) for 
different Q, values, for the case oj^ 1 << (3. The val- 
ues of s(t) increase from zero, corresponding to a pure 
state, to 1/2, corresponding to a completely mixed state, 
with superimposed oscillations. For Ohmic dissipation, 
three time regimes can be identified |30|. The short- 
est time scale is set by t < O (1/oj c ). The quantum- 
fluctuation dominated regime corresponds to O (l/w c ) < 
t < O (1/kT). The thermal-fluctuation dominated regime 
is t > O (1/kT). Our short time approximation yields 
reasonable results in the first two regimes. For t > 
0(l/kT) it cannot correctly reproduce the process of 
thermalization. Instead, it predicts approach to the maxi- 
mally mixed state. 

Figure 2 corresponds to the parameter values typical for 
low temperatures and appropriate for quantum computing 
applications, u> c (3 = 10 3 , with fi = 1.5 • 10 -7 chosen to 
represent weak enough coupling to the bath to have the de- 
coherence measure reach the threshold for fault-tolerance, 
of order 1CP 6 , for "gate" times well exceeding 1/uj c , here 
for oj c t over 10. The leading-order quadratic expansion 
in powers of the time variable t, the validity of which is 
limited to t < O (1/oj c ), is also shown for comparison. 



Figure 1. Schematic behavior of s (t) for different values 
of f2, decreasing from i to iv. 




Figure 2. Comparison between the 0(t 2 ) expansion, i, 
and the short-time approximation, ii. 
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